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I. INTRODUCTION 


In feedback control systems, the plant control is often 
a function of the current plant state vector. A knowledge 
of the state vector is therefore required for the controller. 
In most systems, however, the entire state vector is not 
available through direct measurement. It is therefore 
necessary to estimate the unknown states. An observer is a 
device which can be used to reconstruct all or part of the 
state vector [1]. 

Traditionally, observers and controllers have been 
designed separately. Time and frequency domain specifica- 
tions are satisfied through the selection of controller 
feedback gains as if all of the plant states are available. 
The observer is designed to estimate the plant states by 
proper selection of the observer gain matrix. The observer 
output is used to supply estimates of the missing states to 
the controller [2]. 

In this thesis, after presenting the necessary background 
material, a method for designing observers which minimizes 
a quadratic performance measure is presented. The method is 
also extended to the simultaneous design of the observer and 
PicwconLrollcro ror Simplrerty, only linear, time-invariant 
systems are considered. The design method relies on the use 
of a digital computer algorithm, however, the computations 
are performed off-line. Numerical examples are presented 


in order to illustrate the new technique. 





II. MATHEMATICAL BACKGROUND 


A. TRADITIONAL DESIGN TECHNIQUES 


a adler, linear, 


The dynamic characteristics of. an n 
time-invariant, free system can be described by the linear, 


first order differential equations 


s(t) = Hs(t) (1) 


where s(t) is an nxl vector and H is a constant nxn matrix. 
The response of the system to any non-zero initial condition 
is determined by the elements of the matrix H which are, in 
turn, a function of the system configuration and physical 
constants. By adjusting system parameters, a suitable H ma- 
trix is realized which causes the system to meet certain 
system performance specifications. 

Traditionally, performance is specified as a consistent 
set of time and frequency domain criteria. Typical speci- 
fications are given in terms of rise time, settling time, 
percent overshoot, steady-state error, gain margin, phase 
margin, and bandwidth. These criteria are then translated 
into a control strategy through the use of some relatively 
Simple mathematical tools; generally by trial-and-error 
techniques. Appropriate assumptions are usually made to 
attempt to approximate complex, higher-order systems by 
Simpler, lower-order models that can be handled with greater 
ease. 

This approach, although empirical in nature, and rely- 


ing heavily on the ingenuity and system know-how of the 





designer, is relatively straightforward for simple, single- 
input, single-output systems. Many modern systems are 
highly complex, however, with multiple inputs and multiple 
outputs. These systems do not readily lend themselves to 
Ree acneneeucal techniques alluded to earlier and often 
performance is not easily specified in terms of the tradi- 
tional criteria. It is therefore desirable to have an al- 


ternative method for prescribing desired system performance. 


B. THE QUADRATIC PERFORMANCE MEASURE 
One such index of system performance is the quadratic 


performance measure 


©0 


J = J t*s(t)Qs(t)dt (2) 


O 
which assigns a positive scalar value to the infinite time 
integral of the time-multiplied quadratic form in s(t) and 
QO; the elements of the nxn matrix Q are specified to deter- 
mine the relative contributions to J of the elements of the 
MeCtOm S(t). = Mne value ot the tCime-multiplhier exponent, r, 
establishes the relative contributions to J of the earlier 
aidweacer Vvdtues oF the sclements of s{t). Interpreting J 
as a cost associated with a particular system configuration, 
it becomes the task of the designer to find values for the 
elements of H which are at his disposal that valet the lowest 
Veet or J mise chovce ton i 1s said to be optimal since 
Mea minimizes the eost . 

Examples of systems to which a performance measure of 


the form of (2) can be applied are the observer error state 





equations 
e(t) = Fe(t) (3) 


and the state-error differential equations of a combined 


plant-observer/controller 


=--- ee ee aa apne ° (4) 
e(t) 0 | F} }e(t) 
Equations (3) and (4) will be developed in detail in Sections 


III and IV, respectively. 


C. EVALUATING THE QUADRATIC PERFORMANCE MEASURE 
In Ref. 3 it is shown that if (1) is a stable, linear 
system and if Q is any positive definite matrix, the per- 


formance measure (2) can be expressed as 


J = (-1)7* ris‘ (0)V_,,5(0). (5) 


The real symmetric matrix ee 1s the unique solution of 


the simultaneous linear, algebraic matrix equations 


; (6) 
HV Sele ieesae Oo 
It can be shown [4] that if H is a stability matrix and if 
ORs pOSttive detinite then the matrix Vy 1s also positive 
definite. The argument can easily be extended to show that 


(-1)7**rtv_ is also positive definite. 


the Sseemeclieicestaple if His a stability matrix; 
that is if, and only if, all cigenvalucs of H have negative 
peal parts. 


10 





The details of the derivation of (5) and (6) are given 
in Appendix A. The solution to (6) can be obtained by 
iteratively solving the r+1l sets of n(n+1)/2 simultaneous 
linear equations- or by a technique suggested in Ref. 3. 
The details of this technique are given in Appendix B. 

Evaluation of the performance measure (2) using equa- 
tions (5) and (6) eliminates the need for determining the 
explicit time response of the system but requires knowledge 
of the initial condition, s(0). Since a prtort knowledge 
of this vector is not generally available, the question 
arises as to how to pick s(0) in order to arrive at a mean- 
ingful value for J. 

One technique, suggested in Ref. 5, is to determine the 
average value of J over a linearly independent set of ini- 
tial condition vectors. This is equivalent to assuming 
that the initial condition vector is uniformly distributed 
on the surface of a unit hypersphere and results in a de- 
Sign which is optimal in an average sense. Another tech- 
Nique, suggested in Ref. 6, is to determine the largest 
possible J for every s(0) within or on a hypersphere of 
radius p and select the configuration which yields the 
smallest of these maximum values. This results in a design 
which is optimal in a worst case sense. For the latter 
scheme, it is shown in Appendix C that for a given design 
parameter set {P}, the largest value for J is given by J, 


where 


4since Pie wate rlces V4 are cach symmetric, 


only n(nt1)/2 elements of Leach’ “eed tS} be solved for. 


il 





3,(P}) = CCD rv Pd) (7) 


and where A, [°] represents the largest eigenvalue of [+]. 
Here, the functional dependence of Jy upon the design 


parameter set {P} is stressed. 


D. MINIMIZING THE QUADRATIC PERFORMANCE MEASURE 

To this point we have assumed that the design parameter 
set {P} is unbounded. In general, however, we must allow 
for the possibility of a constrained {P}. There may be 
explicit as well as implicit constraints. The explicit con- 
Straints are those which are the result of the physical 
limitations on the system components while the Imp 1 care 
constraints are those which result from restrictions placed 
on system performance by the designer 

In either case, the constraining relationships must be 
included in the problem formulation. One means of accon- 
plishing this is to add a large penalty to the cost of a 
given configuration if that configuration requires the vio- 
lation of any of the constraints. With this in mind we de- 
fine the optimal parameter set {P"} as that parameter set 
Which minimizes J, ({P}) and the minimum value of Tae 3 
as i Glitters 


* 


J 


i 


3,({P°}) = min J, ({P}), 
. {P} 


Therefore” 


EE ee 
ee ; ; ; 
Since rt! is a constant in any given design specification, 


it docs not contribute to the optimization of the parameter 
set and can therefore be deleted from the performance measure. 


2 


—— ‘es 





Jos minty (OD) Ve CPI, (8) 


This indicates that the following iterative procedure 
can be used to find (P"} once r and Q have been specified. 

1. Select an initial estimate for {P} insuring that 
any constraints imposed are satisfied and that H is a sta- 
bility matrix. 

2. Using equation (6) or the method described in Ap- 


pendix B, find 


jeer lk 
5. Evaluate 


+1 
Cte ae had) | eee Ce) 


4. Through some suitable minimizing search routine” 
select a new set of parameters {P}. If the constraints are 
not satisfied, add a penalty to the cost and select a new 
parameter set. 

5. Repeat steps two through four until some criterion 


for ending the search is satisfied. 


tor the examples wsed in this thesis, the subroutine 
“DIRECT” (NPS Computer Facility) that performs the pattern 
search method of Hooke and Jeeves was used. 


eS 


— 





IIIT. OBSERVER DESIGN 


In the design of feedback control systems it is often 
desirable to have information about plant states which is 
not normally available through direct measurement. It is 
sometimes possible to construct a system called an observer 
which can estimate a linear transformation of-the plant 
state and, if the transformation is invertible, it is then 
possible to reconstruct the plant state vector itself. 

Traditionally, observer design has been accomplished 
through the use of techniques such as those described in 
Section IIA. In this section, techniques for the design of 
full-order and reduced-order observers by using the quad- 
ratic performance measure as the design specification are 


presented. 


A. FULL-ORDER OBSERVER DESIGN 


We consider the time-invariant, linear system 


x(t) = Ax(t) + Bu(t) 


(9) 


Gel eS Gas) 


Wichee VCC rs tne met Statesvector, y(t) 1s the mxi output 
VeCUC@ Ct) ised. checoMtnol Vector, A is an nxn constant 
matrix, B is an nx% constant matrix, and C is an mxn con- 
Stant matrix. We consider also the full-order identity 
observer shown in Fig. 1. The observer state equations as 


developed in Appendix D are given by 
z(t) = (A-GC)z(t) + Gy(t) + Bu(t) (10) 


14 











Pocue ler ulilOncenr Observer. 


LS 





where z(t) is the nxl observer state vector, G is the nxm 
Gbseérver gain matrix (in this case the elements of G con- 
stitute the design parameter set {P}), and where it is as- 
sumed that (C,A) is observable. 


Defining the error state as 


e(t) = z(t) = x(2) 


and differentiating with respect to time yields 
3) 2 2G) = <G2- 
Substituting (9) and (10) for x(t) and Zt Ne respectively, 


and solving for e(t) yields 


oO) STA Ce) el), 


We define the matrix F as 


B= (A-GC) art) 


which yields 


e(t) = Fe(t). Gz) 


We define the quadratic performance measure as 


co 


/ tTe! (t)Qe(t)dt. (13) 


oO 


A 


The system (12) and the performance measure (13) have the 
same form as (1) and (2), respectively. Equation (13) can 


therefore be rewritten as 
Pere rt Lt sy 
J = (-1) 6S (O)V_ 7200) 


where, by Using (6), 


16 





ee oot eee 
VEN > tee 7 Seco) 
ae ee 


Furthermore, the algorithm described in Section IID can be 
k 
applied to find G-, the observer gain matrix which is opti- 


mal with respect to the performance measure (13). 


B. REDUCED-ORDER OBSERVER DESIGN 


We again consider the time-invariant, linear system 


given by 
x(t) = Ax(t) + Bu(t) 
(14) 
ae) een GCs) 
We define the vector x(t) as 
z (t) 
sce & [ESD 
w(t) 
where w(t) is an (n-m)xl vector given by 
w(t) = Dx(t) (15) 


and where D is an (n-m)xn matrix, arbitrarily selected so 


that 
-1 


exists. By proper selection of the state variables C will 


have the form 
Oo eS we 
Penotaerond!)y cConventent chorce tor )) 1s 


D=- [0:11], 
17 








We can find w(t), the estimate for w(t) by constructing the 
reduced-order observer as shown in Fig. 2. For the above 
choice for C and D, the observer state equations developed 


in Appendix D can be written as 


z(t) = (Ay2-GA),)2(t) + [A,)-GAy, + (Ap-GAy»)GIy(t) 
+ (B,-GB,)u(t) (16) 
where 
G5) Sls) caidas) (17) 


and where the matrices A and B are partitioned as 


ay eis oe 5 et 
— 3 — 
Ao, | Ago Dy 
Here, Adi is an mxm matrix, Ay is an mx(n-m) matrix, A543 
is an (n-m)xm matrix, A> 1s an (n-m)x(n-m) matrix, By 1S 


an mx matrix, and B, is an (n-m)xX& matrix. Therefore, 
x(t), the estimate for x(t) (which is, for the case 
C={1! 0] andD=[0: I }, the estimate for x(t)), 


becomes 


ZC ee Gy, (it) 


Pcie -tnewerronr ti) tie eStamate w(t) as 


e(t) © w(t) - w(t) 


HiceoUvStreumimen( 1S )ieand (17) for w(t) and w(t), respec- 


tively, yields 


eine oc) (GS-1) xe). 
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No GN A 


Any Aq to 


u(t) 





PEcurencaes wReaucea-Order Obsérver . 


1S 





Differentiating with respect to time, substituting (14) and 
ClG je Eor Ge) and Zee); respectively, and solving for ecm) 


yields 
e(t) = (> -GA,.)z(*) + 1A,,-GA,, * (A,,-GA,,)Gly(*) 
+ (B,-GB,)u(t) + (GC-D)(Ax(t) + Bu(t)). 


We note that forC = [1:0] andD=[0:1] 


—~ 


GCAx(t) = [GA,, : GA, ,]x(t) 


~DAx(t) = [-A,q -A,>]x(t) 
GCBu(t) = GB,u(t) 
-DBu(t) = -Bou(t). 


MeenotealsO.tiaceroT y(t) =) Cx(t) 
Azy(t) = [A,, : O)]x(t) 
~GA,,y(t) Za [-GA,4 5 O)]x(t). 
These relationships, when substituted into the above expres- 
sion for e(t), yield 
=|| 0 : (A,7-GA, 5) ] x(t). 


We note that the last term in this expression can be written 


as 
= [ 0: (A,,-GA, 5) ]x(t) = 5s (A,,-GA, 5) Dx(t). 


Substituting this relationship into the expression for e(t) 


and combining terms yields 


e(t) = (A,.-GA,>) (z(t) + (GC-D)x(t)) 


20 





where 


PCr emo) XCE yy =e (tc). 


Therefore e(t) becomes 


e(t) = (A, .-GA,,)e(t). 


If we define 


= (Ag27GAq7) eS) 


ie?) 
! 


then 


e(t) = Fe(t). | (19) 


Defining the quadratic performance measure as 


Le «] 


J af tTe! (t)Qe(t)dt (20) 


O 
we see that the system (19) and the performance measure (20) 
have the same form as (1) and (2), respectively, and that 


(20) can therefore be written as 
= veal T 
J = (-1) Coc (O)V 4 200). 


By using (6) we have 


il = = 
Retire viene hem al Selanne 


- : 
FeV, a ViF = =O 


OlGeravalinmeneeadloorrthin described in Section IID can be 
* 
used to find G , the observer gain matrix which is optimal 


with respect to the performance measure (20). 


Cee KANE UE OBSERVER DESIGN PROBLEM 
LOmmiiustrate thevapplication of the®algorithm to ob- 


server design we consider the plant whose transfer function 


Zall 





is given hy 


ACW We ea 
S meet 125° * 475 * 60 





The state and output equations for this plant can be written 


as 
0 1 0 | 0 
x(t) = | 0 o 1 xct) + fol] uty 
=605)-47 -I12 " 
J J 
VAC See OO le  ) 
he ce: 
0 1 0 0 
A = 0 0 1 » B= {0 mee [| 1070) |]: 


OO mums et 7 | 

It can be easily shown that (C,A) is completely observ- 
able. A second-order observer with arbitrary eigenvalues 
can therefore be constructed to provide estimates for the 


unmeasured states X, (t) and xX, (t). These estimates, Ot) 


and X(t), are given by 


eee) 


ee) = WG) = 2Ge) & Gyula) 
x 
3 


We require an initial estimate for the observer gain 
Nahi Gewhrehn yrelds a stable observer. The observer F 


matrix is given by (18) as 
He) ig Ne a) 


Here 


Zc 





> 
te 
> 
NN 
lo | 
—_ 
© 
—s 
se 
6) 
W 


Therefore 
“84 1 
EE ac. -12 


where the observer eigenvalues are the eigenvalues of F and 


are given by the solution to 


A? + (12+g,)A + (12g,+g,+47) = 0. 


We arbitrarily select the initial estimate for the ob- 
server eigenvalues to be -7, -2. The characteristic equa- 


tion is therefore 
(A+7) (A+2) = 0 
or 
A* + 9A + 14 = O. 
Equating coefficients of the like powers of A and solving 
for the G matrix yields 


= 
Gy- 


iricmeno1 Ces fOrsG Mise also satisty all design constraints. 
A parameter in observer design which must be limited is the 
observer speed since, as the observer gains are adjusted to 
SalsewtMw@werrors to dle out more rapidly, the obsrver be- 


comes more sensitive to noise. One measure of the observer 


ie 





speed is the sum of its eigenvalues. It can be easily shown 
that this sum is numerically equal to the trace (the sum of 
the diagonal elements) of the F matrix. 

As our constraint on observer speed we therefore set a 
limit, T, to the most negative value that we will allow the 
Sater tone Mat hixarne@ assume... 0@ observe the effect of 
T on the resulting design we will allow T to take on the 
values -9, -12, and -15. We see that for our initial esti- 
mate for G that the trace of F is -9 and that this G will 
therefore satisfy any of the values of T previously stated. 

The design algorithm of Section IID was used to compute 
the optimal observer gain matrix G”, the corresponding value 
of the performance measure 1G )- the observer eigenvalues 
dy and hos and the worst case initial condition error vector 
for various choices of the design specifications Q and r in 
the performance measure (20) and for various choices of the 
constraint T. The plant and observer were then simulated 
on the IBM 360 digital computer for each design specifica- 
tion using the corresponding G and the worst case initial 
condition error. 

Table I shows the effect of varying the design parameter 
r on ce Wee). and the observer eigenvalues hy and ho. For 


these cases the matrix OPiS xed at 


and T is fixed at -12. Note that in each case, the value 
se 

FOr 1 iiPZeGe oi CemeneCast Gace must be greater than or 

equal to T we see that 


24 





iz-- 2, 2 - 12 


is satisfied only for g, less than or equal to zero. The 
fact that g, is always zero indicates that the constraint 
is always active. 

The effect of varying the parameter r on the observer 
error time response is shown in Fig. 3. Note that as the 
value of r is increased, the estimate error for both states 
reaches zero sooner. For the estimate error for the state 
X(t), however, this is at the expense of larger error 
values in the early part of the operating interval. 

In Table II the effect of the design parameter 444 07 


# 
G , J(G@ ), A, and A, is shown. Here, the value of r is 


2 
fixed at zero and the value of T is fixed at -12. The ma- 


trix Q has the form 


Again we note that the constraint T is active. 
The observer error time response is shown in Fig. 4. 
The effect of increasing 4), 0n the plant states is similar 
to that of increasing r, although somewhat less dramatic. 
iimelabbe shi tisthie ctfect of the design parameter do. on 


x a 
cr (Ga). \,, and A, 18 Shown. Again, the value of r is 


7 
hmixea at zero and the value of T is fixed at -12. The ma- 


trix Q has the form 


OF 


© 








x fe 
hab lew temeitcwe =fCeCrsot 1 ON. Gear) (G hy and do. 


0 
Q= |p ‘ , T = -12 






-6.0-j1.4 
-6.0+51.4 


# K 
Table II. The effect of q,, on G, J(G), A, and A,. 


r= 0, T= -12 
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Once again we see that the constraint T is active. 

The observer error time response is shown in Fig. 5S. 
The effect of do2 on the response is quite different from 
the effects of r and 441 aS we might expect. By placing 
greater emphasis on the cost associated with states which 
represent higher derivatives we do not allow them to devi- 
ate as far from zero and thus tend to slow the convergence 
process down. Hence, as d52 increases, the observer error 
states both reach zero at later values of time. 

It should be pointed out that all observer configura- 
tions were simulated for a period of two seconds for the 
sake of uniformity in the graphical presentations. If the 
response of Fig. 5 were continued for a greater length of 
time, the curves would all eventually approach zero. 

In Table IV the effect of the design constraint T on 
G 5) hy and do is shown. Again, ris fixed at zero 
Sic Quls! bixed at 


1 0 
Q = 


By noting the values for g, we see that the constraint T is 
always active. 

The observer error time response is shown in Fig. 6. 
We see that the effect of making T more negative is to im- 
prove the shape of the response at the earlier values of 
time. Note that in the case where T has the value -9 that 
the observer error in state x5 (t) actually gets worse for 


early values of time. 


Zo 











* * 
Table JII. The effect of do2 on Gey (Gy) ST 
r= 0, T = -12 
x * 
htptemi Vim Mieworticce 0% | on G.; ¢{G ), hy and do 


Ses0) = 
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IV. OBSERVER/CONTROLLER DESIGN 


A. TRADITIONAL DESIGN TECHNIQUE 

The traditional method for the design of feedback con- 
trol systems utilizing observer estimates of unmeasured 
states is to specify desired plant performance in the manner 
described in Section IIA as if all states are available 
through direct measurements. A feedback network or controller 
is then constructed which causes the closed-loop system to 
satisfy these performance specifications. The controller 
Bias available as inputs either the observer estimate of the 
entire state vector or the measured plant states and the ob- 
server estimate of the remaining states. The observer is 
designed to be stable and to have a characteristic equation 
with eigenvalues which satisfy the somewhat nebulous require- 
ment that the observer by "slightly faster than the plant" 
i2]; 

In this section, an alternative method is described in 
which the observer and controller are designed simultaneously 
using a performance measure similar in form to (2) in Sec- 
mton LIB. This results in a configuration in which both the 
controller design and the observer design are optimal with 


respect to the specified performance measure. 


B. FULL-ORDER OBSERVER/CONTROLLER DESIGN 
The full-order observer for the linear, time-invariant 


system 
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x(t) = Ax(t) + Bu(t) 


(21) 

DACe = Coste) 

is, as shown in Appendix D, 
z(t) = (A-GC)z(t) + Gy(t) + But). (22) 


We define the control law as 
Wie) = hx) 


where K is a constant £xn matrix. 

It can be shown [2] that the stability of a closed-loop 
system using an observer output as a substitute for the 
plant states in the control law is unaffected by a stable 
observer. The control can therefore be approximated by the 


control estimate, Ge) where 
u(t) = Kz(t). (23) 


The structure for this full-order observer/controller is 
shown in Fig. 7. 
Defining the error state, e(t), as 
A 

Cee) = 21h) > x(t) (24) 
and differentiating with respect to time yields 

e(t) = z(t) - x(t). 
Substituting (21) and (22) for ae) and FAG) respectively, 
and solving for e(t) yields 


e(t) = (A-GC)e(t). (25) 
Substituting (23) into (21) yields 
x(t) = Ax(t) + BKz(t). (26) 
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Figure 7. Full-Order Observer/Controller. 
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Solving (24) for z(t) and substituting into (26) yields 
x(t) = (A+BK)x(t) + BKe(t). (27) 
Defining the augmented state vector, &(t), as 


A | Xt) 
GE ele ea 
e(t) 


and using (25) and (27), it is easily shown that 


A+BK — BK 


E(t) = EC) - (28) 
0 A-GC 


The performance measure can be defined as 


c (29) 
J = [ t™[x'(t)Qx(t) + e' (t)Re(t)]dt. 


O 
This performance measure is chosen because it accommodates 
the specification of convergence properties of both the 
closed-loop system states and the observer states and im- 
plicitly limits the expenditure of control effort. Equa- 
eron (29) can be written in terms of &(t) as 
Qo 


J - | tTe! (t) E(t)dt. (30) 
Q 0 


|x 


Equations (28) and (30) are in the same form as (1) and 
(Z), respectively. It has been shown [2] that the matrix 


¥, defined as 


36 





has as its eigenvalues, the eigenvalues of A+BK and the 
elgenvalues of A-GC. The matrix ¥ is therefore a stability 
matrix if both A+BK and A-GC are stability matrices. It 


can be shown that the matrix 


Q 
OR 


{oO 


1s positive definite if both Q and R are positive definite 


[7]. Therefore (30) can be rewritten as 
2 toa ol eed 
J = (-1) re& (O)V_ 8 (0) (31) 


where, by using (6) 


AtBK BK Y! A*BK BK 
; eee Vegy V4] ie = ae SS Zrocee ls 
A+BK BK 1! A+BK BK Q 0 
7 ee eee | F 
0  A-GC 0 A-GC oR 


Furthermore, the algorithm described in Section JID can 
be applied to find K and Gs the feedback and observer 
gain matrices which are optimal with respect to the per- 
formance measure (29). 

A slight modification to the design algorithm causes a 
Significant reduction in the number of computations re- 
quired on each iteration. To see why, we note that if the 
PUSehVer E€Stilate indtial Condition vector 1s aSsumed to be 


Mic null yeetor, that 1s 


Ot 





z(0) = 0 
tne igmer Tole. 5). ¢ (0) becomes 


SO) =e) 


and the augmented state vector &(0) becomes 
BO) “=| =e>= (32) 


Partitioning the symmetric matrix Ur as 


V 


a 


—12 


Neel = 
T 


we. oy, 


where each element is an nxn matrix, then substituting (32) 


and (33) into (31) yields 
(34) 


r+ 


J= (-1)7 7 etx! (0) 1V,4-Vyo-Vy 9 *¥9 2) x00) - 


‘= 12 —12 


We now define the nxn matrix ved as 


_ ; : i 
Te 2 Wa Gin Se i) 
Therefore 
Eel ah 
J = (-1) eS (O)V 1X00). (36) 
The matrix Ne has nectgenvalues, while the matrix en 


hasecimeigenvalues. “We see, therefore, that the effort 
Required in finding Jy can be significantly reduced by using 


(35) and (36) in step three of the design procedure. 
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C. REDUCED-ORDER OBSERVER/CONTROLLER DESIGN 


For the linear, time-invariant system 


x(t) = Ax(t) + Bu(t) 


(37) 
ne) = Gale) 

we define the control law 
u(t) © Kx(t) (38) 


where K is a constant &xn matrix. We define the vector 


Tes 
< 
pie 
ct 
Need 


x(t) 


where w(t) is an (n-m)x1 vector. We let 
WCC) DC cy (S53 


where D is an (n-m)xn matrix, selected so that 


exists. By proper selection of the state variables, C will 


have the form 
C= { 1:0 ]. 
A notationally convenient choice for D is 
D= [ Of1 ]. 
hor chis cho1eee of Cand D, w(t), Giewestimate for w(t); 
can be obtained by constructing a reduced-order observer, 


as developed in Appendix D, whose state equations are given 


by 


ESS, 





z(t) = (A,,-GA,.)zZ(t) + (A,,-GA,1*(A,2-GA,2) GS) y(t) 


: (40) 
+ (B,-GB,)u(t) 
where 
w(t) = z(t) + Gy(t). (41) 
The matrix K and the vector x(t) can be partitioned as 
m 
x (t) 
K= [Kiko] , x(t) = |=5--- 
— x(t) 


where, for the stated choice for C and D, 


y(t) = x(t), w(t) = x°(t), w(t) = x°(t). 


The control law, given by (38), can therefore be approximated 


by 
A x" (t) 
ice) = K | pees 
or 
u(t) = [K™:K°] ores seas | 
z(t) + Gy(t) 
or 


MON ake te Koz it) > KGy(t). (42) 


The structure for this reduced-order observer/controller is 
shown in Fig. 8. 
Med oLincmtnceGusemyat hon error, €(t), as 
A “A 
CO) a= WCC eat t)., (43) 
Substituting (41) for w(t) aio jeson wit) yields 
CCH oe oC eet Gy it) s- Dx (tt). oy) 
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Eigure 8. 
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Reduced-Order Observer/Controller. 
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Substituting (37) for y(t) and rearranging yields 
CC.) ez (SC- 0) x(t): 
Differentiating with respect to time yields 
e(t) = z(t) + (GC-D)x(t). 
Substituting (37) and (40) for x(t) and AG Besmee eiVewy.. 
and solving for e(t) yields, as shown in Section IIIB, 
2100), = (A,-GA,,)e(t). 
Solving (44) for z(t) + Gy(t) and substituting into (42) 
yields 
a ill O 
MGe aS ey CO) ate Ct) Dxee 
We note, however, that 
Dx(t) = x°(t) 
and 
VAG) een Ce 
Therefore, 


u(t) 


Kx(t) + Kve(t). (46) 
substituting this estimate for u(t) into (37) yields 

x(t) = Ax(t) + B(Kx(t) + K°e(t) 
which, after rearranging, becomes 

x(t) = (A+BK)x(t) + BK°e(t)). (47) 


He detine the avigmented state vector, €(t), as 





Using this definition and the expressions for x(t) and e(t) 


given by (47) and (45), respectively, yields 


A+BK  —BK™ 


E(t) = E(t) . (48) 
Q A 


The performance measure can be defined as 
(49) 


oo 


[> f tT rx! (t)Qx(t) + eT (t)Re(t) }at 


O 
which can be rewritten in terms of the augmented state vec- 
EGG ot), as 


Z Q 
J = f te! (t) E(t)dt. (50) 
0 0 


Le) 


|x 


Equations (48) and (50) are in the same form as (1) and 


(2), respectively. Furthermore, the matrix 


A+BK  BK® 


0  Ayy-GA 


a2 12 
ioeieStaD Eli ty Tathixeis both A+BK and A,.-GA,, are stability 
Matrices. Also, the matrix 
Qo 
0 R 


PP OtiMmenconin ttc mt bounmO. and Rare positive definite. 
The performance measure (50) can therefore be rewritten 
as 


Pee (Ove, 2 (0) (51) 


J= (-1) —r+]— 
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where, by using (6), 


AtBK BK] # AtBK —BK® 
V V - 
Suse S70 i 7 =s 
22 2 oD Gee s0. 
: s = 1,28 
A‘BK  BKo ] ¢ A*BK —BK® Q 
Ss ia 
OO Bg a, Uh so Boy MES, u 


Furtnermore, the design procedure of Section IID can be 
applied to determine K" and ene the feedback and observer 
gain matrices which are optimal with respect to the speci- 
fied performance measure (49). 

Once again, it is possible to reduce the number of 


computations required for each iteration of the design 


algorithm. In this case, we set the observer estimate, 
w(0) TU aiheMeysehdcme (Oj “Dx ( 0) ae ineretore —€(0) be- 
comes 
x (0) 
| es 
-Dx (0) 


(52) 
T 
ee es STE nes tars 
-Dx (0) -Dx (0) 
Pavel elon enew( 2n-m)x(2n-m) symmetric matrix ME se as 
yo. 2 [ott 742 
—r+] T 

ui 22 
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where Vi4 


V59 is an (n-m)x(n-m) matrix and substituting into (52) 


yields, after expanding 


is an nxn matrix, V4 1s an nx(n-m) matrix, and 


(53) 
a teal T ae ie A 
J= (-1)° “rix (O)(V,,-D Vi -V¥y2D*D V5 DIx(0). 
We define the nxn matrix V! as 
—r+1 
A fea ie T 
eee me n> DDD oo 
Therefore, (53) can be rewritten as 
2 jis el T 
J = (-1) eG (O)V x (0). (55) 
The matrix ta has n eigenvalues, while the matrix ae 


has 2n-m eigenvalues. We see, therefore that the effort 
Peqiite dn inemade te ria name Jy in step three of the design pro- 


cedure can be reduced by using (54) and (55). 


Doe EXANPEE OBSERVER/CONTROLLER DESIGN PROBLEM 
To illustrate the application of the design algorithm 
to observer/controller design we consider the plant whose 


transfer function is given by 


CS) 1 


U(9) Se last 20 


The state and output equations for this plant can be written 


as 
; 0 1 0 
(UE) >= eaGt) Ge) 
-5O -15 
Meese Ox Ct). 
me le: 
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A = B= wey owe oe, [ead tea 0 a 
-50 -15 1 
Then 
ike ks] 
0 1 to 
A+BK = 7 
-50 -1S5 1 
0 i} 
k,-50 k,-15 
Also 
Fr As2 ~ GAr2 
where 
Mog @ (Lleol) aCe aligl 
Therefore 
a eee ale 


NewmmMust "seleet an Phixi at. decimate for both K and G 
WibtelincavsccetNeunadtnices AtBK and Foto be stability matrices. 
As in the case of the observer alone we must constrain sys- 
tem speed. In this case, however, we sum the eigenvalues 
of the closed loop system. It is easily shown that this 


SUIS eIVven Dy the trace of the matrix 


AtBK BK” 

: Se, 
which is given by 

kK, eae) a su 





We arbitrarily select the trace constraint T to have a value 


of -25. Therefore 


= Zoe k, ee a 


We see that the open-loop system is stable. We therefore 
choose an initial estimate for K of [0 0]. We must now 


choose a value for g which satisfies both 


Coes wom 2 
and 


mee oO aac 0 
Or, combining the two inequalities, 
SL Dee eae ne 
NeEmenoOse tO OULrscStLIMmate Of GG, thererore, 
Gee ele Onl 


The design algorithm of Section JID was used to compute 
the optimal controller feedback gains ‘a and observer gains 
cu the corresponding value for the performance measure 
FOR eee the plant and observer eigenvalues his d. and hes 
and the worst case plant and observer estimate error ini- 
PUacCONalelolns L£On Various choices of Q and r in the per- 


fOrmanece measure (49). The Constraint T was fixed at -25 


and the design parameter R was fixed at 
ee th) 
inpall cases. wine plant and observer were simulated on the 
JBM 360 digital computer for cach pé@rformance measure spe- 
* * 
Gvived Using Koy G 5; and the corresponding plant and esti- 
Naceeworst Case initial conditions. 
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Table V shows the effect of varying the design parameter 


* * * * 
Pete eso (Kk 5G ), his r el ha JU te ole Ve RIES Melnkos (0) 


yd 5) 


matrix is also fixed at 


Figure 9 shows the effect of varying the design parameter r 
on the closed-loop system response to the worst case initial 
conditions. We see that increasing r improves the system 
settling time in both states at the expense of larger values 
of x, (t) iene sank warts oOrstnenOperatgno si mtcryal. 

Table VI shows the effect of varying the design parameter 
4,1 on K, Co Oe 6 4), hy: hos and Az. In this case r is 
fixed at zero and Q has the form 


Siew 


Q = 
0 1 


Figure 10 shows the effect of varying 3,4 0n thevelosea-) oop 
system response to the worst case initial conditions. We 
SeCCetladt. ds the value ot qi ommMenheasecdmmtene Settling time 
improves for both states at the expense of larger values of 
x, (t) in the early part of the operating interval. 

Table VII shows the effect of varying the design parame- 
ter qj, on K", e, J(K’,G’), Mos hos and Az. i ul sacase 
Eebouirxed At zero and Q has the form 
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ss = 





apie V. 





= S=4 


-11.5+j11.6 





Table VI. 


The effect of q 


and Az ewes ee 25) oR = |i). 


ee 
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Figure 11 shows the effect of varying d>> on the closed-loop 
system response to the worst case initial conditions. We 
see that as 459 is increased, the system becomes more slug- 


gish, requiring longer settling times, as we would expect. 
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V. CONCLUSIONS 


A new technique for observer and observer/controller 
design which requires the minimization of an integral of 
time-multiplied-square-error (ITSE) performance measure has 
been developed. An algorithm which is well suited to digi- 
tal computation and which is especially useful for higher- 
order systems has been presented. An example of a computer 
program which implements the algorithm for the case of an 
observer/controller is shown in Appendix E. The computa- 
tions are performed off-line with a typical computation time 
Sree miyey seconds. 

It has been demonstrated through the use of examples 
that modification of the time response of an observer or 
aeclosed-loop system wLilizing an observer/controller to 
initial condition error can be achieved by a suitable change 
in the appropriate design parameters of the ITSE performance 
measure. Furthermore, it has been demonstrated that the 
qualitative change in the system time response reflects 
diGeculysthe Gudlitative Change in the performance measure 
design parametcr. 

Tt is recognized that for sccond-order and, perhaps, 
third-order single-input, single-output systems the rela- 
tionships between traditional design specifications and 
System transtér functions are simple and relatively easy to 
Hise In tmese cases, the design procedure suggested above 


1s generally not efficent. However, when dealing with 


=e 





higher-order, multiple-input, multiple-output systems, the ad- 
vantages of the new technique become more significant. For 
these systems, the design specification and system transfer 
function relationships are highly complex and the traditional 
methods are therefore extremely difficult to use. However, 
the new technique offers a relatively straightforward ap- 
proach to design which does not become more complex as the 
system order, the number of inputs, or the number of outputs 


increases. 
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APPENDIX A 


EVALUATION OF THE INFINITE TIME INTEGRAL f{ t*sl(t)Qs(t)dt 
O 


For the system of linear, first-order differential 


equations 


s(t) = Hs(t) (1A) 
if there exists a scalar function 
V,(s(t)) = s' (t)V,s(t) (2A) 


subject to the conditions 


1) V,(0) = 0 


Zz) V,(s(t}) 20g, EOrman s(t) 3 


ge, 


3) Vy(s(t))> ©, as [|s(t)|| > © 
AyeV, (sile)) <0, for alll s(t) + 


5) Vy(s(t)) = 0, s(t) = 0 


io 


then the system (1A) is asymptotically stable [4]. 
Since V,(s(t)) 1s a scalar, it 1S equal to its transpose. 
Hence, from (2A), we see that the matrix Vv, Ps esymmetric, 


We now consider 


V(s(t)) = 5) (t)Vys(t) +s’ (t)V,8(t). 


SuUpStrtutane (1A) for s(t) yields 
V,(s(t)) = (tiv, s(t) + s"(t)V,Hs (t) 


which, after rearranging terms becomes 


7 


V,(s(t)) = s'(t)(H'v, + VjH]s(t). 


We now define the matrix Q as 


qo 5 


and note that since V, is symmetric, then Q must also be 


iilcer’ ec... liemet one 
VGC) = sstGyeste). (3A) 


At this point we note that equation (2A) insures that 
conditions one and three are satisfied. Condition two is 
satisfied if V, is positive definite and conditions four 
Micmeivcmarcos@bistiled tf ONS positive detinite. 

Integrating the left hand side of (3A)} over the finite 


Pimearnucmyal.. | eweyieclds 


iT T 
[¥yecoae = f avy (sce) 
O O | 


it 
[¥,cscep ae V,(s(T)) - Vy(s00)). 
0 


If (1A) is asymptotically stable, then 


U7. 


Lim V,(s(T)) 
T>0 


Therefore, 


oo 


J ¥ycscen pat 


O 


0) 





Or, from (ee 


co 


J s'(t)Qs(t)dt=s (0)V,s(0) | (4A) 
O 
where Vi is the solution to the algebraic matrix equation 
Bev, = ll = oe (5A) 
BO iL 


We now consider 


aglts(t)V,s(t)] = s'(t)V,s(t) * 18) (t)Vs(t) + s°(t)V,8(t)] 
= V,(s(t)) + tV,(s(t)). 
Hence 
J ditts™yyysct] = if V,(s(t))dt + | tV, (s(t))dt. 
O O O 


Mepean Ocesnown that wit s(tj ers the state vector of an 


asymptotically stable, free system, then for a finite n, 


- d{t"s (t)V,s(t)] = t™s (t)V,s(t) =D 
O O 
Therefore 
0 = if V,(s(t))dt | tV, (s(t))dt 
O O 
or 
if tV, (s(t))dt ~ | V, (s(t) jdt 
O O 
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or, substituting (3A) for V,(s(t)), 


Co fee) 


f tsTcegsceyat = - if V, (s(t) )dt. (6A) 


O 0 


Now suppose there exists a scalar V,(s(t)) such that 


EGC) 2 = iso (7A) 
and 

V(s(t)) = V,(s(t)). (8A) 
ier 

_oT T 

Vo (s(t)) = s (t)Vos(t) + s (t)V,s(t) 
Or 


V(s(t)) = s (t)(H'V, + VHis(t). 
However, from (8A) and (2A) we have 


V,(s(t)) = s° (t)¥,s(t). 
Therefore, we see that 


HV, + VoH = V 


Z —Z— 1 


where, since V, is symmetric, then V, is symmetric. From 


i 
(8A) we have 


oo 


- f vycsceyyae 


6) 


- f Wescenat 
0 


oD co 


- if V, (s(t ))dt = J dv,(s(t)) 


0 oO 
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co 


-f V, (s(t) )dt 


0 


V(s(0)) 


co 


-{ V, (s(t) )dt 


.@) 


s'(0)V,s (0). 


Substituting (6A) yields 


co 


-f ts (t)Qs(t)dt = s (0)V,s(0) (9A) 


0 


where V, is the solution to the algebraic matrix equations 


HV, + VoH = V, 
(10A) 


HV + Vii ss 


The argument can be continued for higher powers of t. 
Thus for a stable system (1A) and a symmetric, positive 


derimite ©. 


co 


vee b _ lear dh T 
if tes aCe s Cedar = (-1) eS: (O)V_4 1800) (11A) 
O 
where a is the solution to the simultaneous algebraic 


Matrix equations 
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APPENDIX B 


A METHOD FOR THE DIRECT SOLUTION OF THE SYMMETRIC MATRIX ere 


A procedure for generating Viet from the system of 


Simultaneous linear algebraic matrix equations 


it - 
2) ier eee 


; : 
HV, + Vb = - Q 


es described in Peyeomenice 3. The procedure, which is es- 
sentially repeated here with minor amplification and nota- 
tional changes involves the dimeee sOlurton for Vea and 
bypasses the need for solving the (r+l1)(n+1)n/2 simultaneous 
equations. 
We first form the matrix S as follows: 
1) Form the symmetric array of integers 
ee ey, 
dy Oe poe 


i z 4569 
r= pi a) alee 


“*n(nt1) /2 


Where Ti is the dimension of H, Q, and ere 
2) Form the n intermediate square matrices, Si; 
tH152,...,n wiose cicments are equal to one and whose order 


mS ens | if 2. 


62. 





3) Proceed along the first row of I.. and for each 


integer, j, that is missing from the sequence 
tieza... intly/2), set all entries in the yth row and column 
of S4 equal to zero. 

4) Proceed row by row through the H and Sy matrices 
Simultaneously. For the first non-zero entry in the first 
row of S4> enter the value hy: for the second non-zero 
entry in the first row of Sy> enter hyo» and so forth until 
Pimeomerwes Ot i itave been tised Co replace all the non-zero 
entries of S,. 

5) Multiply the first column of Sy by two. 

6) Proceed along the second row of ix as in step three, 
making the appropriate entries in the S, Matrix, 

7) Repeat step four for the S, Ta (ede, 

8) Repeat step five for the second column of the S5 
Wa td X:, 

9) Continue the procedure until all of the S; matrices 


have been generated. 


10) Form the S matrix from 


n T 
S im 421 55 
We next form the vector Lose | from 
Le> ee) 
ell ss S q 


where eee) and q are the leading diagonal elements of the 


Symmetric v 1 MO Mite rtees. nespeceively. That 15), 


T+ 
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Wel 





Inn 


The matrix Vee can now be generated from the vector Were 


in a simple manner. 
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APPENDIX C 
a emia ice 
EVALUATION OF THE QUADRATIC FORM s (0) (-1) reV4 500) 


Let x be an n dimensional, unit length vector. That is 


x! x = ie LG) 


The n dimensional vector y whose direction is the same as 
enateor x but whose length, p, 1s arbitrary is therefore 


given by 
Y= PX - (2C} 


We consider the quadratic form 


c= yay (3C) 
where A is symmetric. Substituting (2C) for y yields 
Gas px Ax. (4C) 


etvus Cons train y to lie on or inside an n dimensional 


hypersphere of radius R. That is, 


O<p <R. (SC) 


We assume we have no further a prtort knowledge of the 
feCuGte, sancatherefore wish to consider its worst case value, 
that is the value of y which causes C to take on its maximum 
value for a given A. 

Miewsec: trom (4¢C) that € is maximized for a given A when 
6 1S a maximum and when the scalar quantity x Ax assumes 1ts 
maximum value. From (5C) we see that the maximum value for 
o 1s R. Therefore, 
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Jt 
Serpe 
(Ces max{x Ax}. (6C) 
We define the quantity, q as 


q 2 max{x! Ax} (7C) 


It can be shown [8] that the matrix A in the form 


x" AX (8C) 
where x! x = 1 and al = A, has n real eigenvalues, d. and n 
corresponding real eigenvectors Xs ii. 2 ess Furkhes= 


more, the maximum value of the form (8C) is the largest 


eigenvalue of A. That is 
Ax} = A fA). (9C) 


Also, the vector which gives this maximum value is the eigen- 
MeCcbor COrresponding to the largest eigenvalue. 


Substituting (9C) into (6C) yields 


= 2 
eee alts veers He (10C) 


Thus, C is the value for C in (3C) which represents 


max 
McwwOuseecasc Value tor y Lor a given A. We wish to mini- 
mize C as a function of A for this worst case condition. 
Pinece nals a itxXedsduantity, 16 plays no role in the mini- 
mization of C and we can therefore assume that it is equal 
Boone tilsedSeCquivalenesto assuming that y lies on the 


unit hypersphere. 


We are concerned with minimizing the quantity 


J = s'(0)(-1)7'riv.,,s(0) (11C) 


eo 
—_——_ 
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rei 
Yr 


where the matrix (-1) eg is symmetric and positive 


*—rt+l 


definite and where s(0) represents the worst case initial 

COncdittonwot thewstate veceor S(t}. We can assume, without 
ligeseot generality, that S(0)} Ties on the unit hypersphere, 
From the argument above, we see that Jas the maximum value 


for J for a given parameter set {P}, is given by 


i. 2 en en 


1 max *—rtl (12C) 


Fide na te tne Worst case inttval condition vector, s{(0), is 


Pete San 


the eigenvector which corresponds to xX Weep 


max 
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APPENDIX D 
OBSERVER THEORY 
oe bois NEED FOR OBSERVERS 


Given a free plant characterized by the state and output 


equations 


x(t) 


y (t) 


Ax(t) ; (1D) 


Gxt) (2D) 


Mmicmcex(bjels the nxl state vector, y(t) is the mxl measure- 
ment vector, A is a known nxn matrix, and C is a known mxn 
matrix, we wish to know the value of the state x(t). If C 
were an nXn matrix representing n linearly independent mea- 


surements, then Gea would exist and x(t) could be found from 
= Ae 
Gees Ge a) (CC) (3D) 


This is not usually the case, however, so we will assume 
that there are fewer linearly independent measurements than 
there are states. That is, we will assume that m is less 
ena MN. 

Since Cm IGe=enOENeors te LOn Wie le ss than n, an alterna- 
Min CmcciciCmMnUls t Oc LOMnd sEomdetermine X(t). One method of 
ES UlNg omen tNnewstate vector 15 tO GoMstruct a model of the 
Pidiirand. operate 1t Under the same conditions as the plant. 
Such a model is called an observer. The following sections 
a@velop the staté Equations for two different types of 


observers. 
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B. FULL-ORDER OBSERVERS 


Let the observer state equations be given by 
2(t) = Az(t). (4D) 
The solution to (1D) is given by 
x(t) = e4'x(0) (SD) 
and the solution to (4D) is given by 
z(t) = e&tz(0). (6D) 


Riemer Orc ei 70) = x(0) 3) the Gbserver State will be iden- 
Phere tO tne plant State for ali t greater than or equal to 
oliee. hOWCVCh, 1m peneral z{0) doeésenot equal x(0), 
the identical dynamics of (1D) and (4D) would result in the 
persistence of the initial condition error and the observer 
state would not provide an accurate estimate of the plant 
SEace . 

An alternative observer is one which contains not only 
the plant model but also a feedback forcing function made 
up from an error signal. Since the only accurate indica- 
MICievcmtavemor the States 15 the measurement vector y(t), 
fetolceassumestnat tne ¢rrer Sipnal represents the difference 
between the plant measurement vector and an identical 
Sieastnemen: “Or the Opserver state vector. Specifically 
we assume that the model and the feedback error signal are 


represented by the observer state equations 


Ate Az) eGiy (tt) s6Cz(t)) (7D) 


dmé {feedback is provided by the’ nxm observer gain matrix G. 


Levieeine (7B) in terms of x(t) and z(t) we see that 
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Game Ch) Gecx(t) - 2(t)). (8D) 


Pectin ewe eCLroLr Stateveur) as 


e(t) £ z(t) - x(t) (9D) 


and substituting into (8D) yields 
z(t) = Az(t) - GCe(t). (10D) 
It is clear, therefore, that by using the proper choice for 
Eomecnecscunon, C(t), can De Wsed to drive Z(t) toward x(t). 
More precisely, differentiating (9D) with respect to 
time yields 
elGe)  BGe = ss(C)- 


Substituting (1D) and (10D) for x(t) and Ze respectively, 


yields 

eG) = hwGy = Gea > esis) 

el) = Male) = sco eae) 

e(t) = (A-GC)e(t). (11D) 
That is 

a(t) = Fe(t) (12D) 


where we define F as 


2 


FEA - GC, (13D) 


The solution to (12D) is given by 


e(t) = eXte(o). (14D) 


memscchmincrerOre, that 1f theve,ecnvalues of I all have 


De Cite VO we Be panes, stmcm for ad finite ¢(0)., 
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lim e(t) = 0. 


t+ 
It is shown in Ref. 4 that the eigenvalues of the ob- 
Selver Calpe posteioned arbitrarily 1£ the pair (C,A) is 


completely observable. That is if the matrix 


has rank n. From (13D) we see that the eigenvalues of F are 
controlled by the selection of the observer gain matrix G. 

The above argument is easily extended to a plant which is 
Daivcimuyacrie Input ur) ines plant State anc outpuc, equa 
tions become 


x(t) = Ax(t) + Bu(t) (15D) 


ea Coa ey) (16D) 


and the observer state equations become 
ZG ec e eeiCS) wee zCC) ey Bute) (170) 


where u(t) is an £xl control vector and B is a known nx 


Matrix. 


fo REDUCED - ORDER OBSERVERS 

The full-order observer possesses a certain amount of 
redundancy since we estimate some states which are essentially 
being measured directly. We therefore desire to construct 
an observer which estimates only those (n-m) states which 
are not being méasured diréctly. We start with the state 
and OUtplt equations 


x(t) = Ax(t) + Bu(t) (19D) 


fell 


Oo 





y(t) = Gx(t). (20D) 


We define the new state vector Ge) SOmtaat tS first m 
components are the measurement vector y(t) and the remain- 


ing (n-m) components are the vector w(t). That is 


tc 
x(t) 4 xt) (21D) 


w(t) 


If w(t) is selected so it is some linear transformation of 


the state x(t), namely 
w(t) = Dx(t). (22D) 
Then, in view of (20D) and (22D), (21D) can be written as 


€ 
Gti ee eC). (23D) 
D 


Furthermore, if D is chosen such that 


=) 


exists, where 


ENen. tne State x(t) can be found from rele) - Namely 
eaten 
ese)! Sa exile 


Assuming that M} exists, the task becomes one of obtain- 
Hie ealiweotanatesror w(t). eWritime the state equations for 


Baot ) we have 


Vz 





Cx(t) = CAx(t) + CBu(t) 


Se 
rc 
ct 
<7 
N 


: : (26D) 
w(t) = Dx(t) = DAx(t) + DBu(t), 
inate is 
x(t) = MAx(t) + MBu(t), (27D) 
SubSstatucinge (25D) fo, x(t) yields 
X(t) = MAM 1X(t) + MBu(t). (28D) 
Defining 
AYE Mam’ BY = MB 
then 2 Reesor Ket) (29D) 





Ae AS 1s t B 
me) = sale | ni = anes + ee u(t) 
w(t) Any : Ad w(t) B> 
(30D) 
Therefore 
Cae Crate) ees wdtect BY u(t) (31D) 
and 
Ce aa ct) Byu(t)y. |) (52D) 


DInGe Viitjeis a measured quantity, 1fS time derivative, 
Ceca Cemiectehniined=:  biewInpuc u(t) Gan also be con- 
SLdened as sa (known quantity. iieteroreemsolvine (SLD) for 


Aj W(t) yields 
' = : 2 ae 7 ' 
Ajyow(t) = y(t) - Aj y(t) - Biult) 
where the é¢lements on the right represent "known" quantities. 


ES 





The quantity Ay w(t) can therefore be used as the measurement 
vector for a full-order observer for the plant whose state 
equations are given by (32D). We therefore write the ob- 
server state equations as 

w(t) = W(t) BMS yalGe) ae JS (Ge) 


(33D) 
+ G(y(t)-Aljy(t) -Bru(t)-At,@(t)) 


where G is an (n-m)xm matrix, w(t) TS thewestimate Or Willi, 
and where the first three terms represent the plant model 
and the last term represents the feedback term. Rewriting 
(55D) yrelds 
w(t) = (A},-GA},)w(t) + (A3,-GA},)y(t) 
(34D) 
+ (BS-GB)u(t) + Gy(t), 

It has been shown [4] that if (C,A) is observable then 


(A! 32) 1s also observable. Therefore G, the observer 


pe 
gain matrix, can be selected so that (A,5-GA1 >) has arbitrary, 
eigenvalues. 

Since differentiation is usually troublesome to perform, 
it is desirable to eliminate the term y(t) EO 054 Di lO 
accomplish this we subtract the term Gre) from both sides 


Sa loa) to get 
w(t}-Gy(t) = (A3,-GA!,)w(t) + (A$,-GA!,)y(t) 


+ (B}-GB!)u(t). (35D) 


Adding and subtracting (A5>-GA) 5) Gy (t) from the right side 


or 35D) yields 
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w(t)-Gy(t) = (A},-GAt,) (w(t) -Gy(t)) 


cee, a7 SA) SY) 


+(BS-GB!)u(t). (36D) 


We now define 


z(t) 2 w(t) - Gy(t). 


Equation (36D) can therefore be rewritten as 
AGES (Cai OY lO} oP Gay Ei Civ eins Bu [02s 
+ (B,-GBj)u(t) (37D) 


where w(t) can be computed from 


w(t) = z(t) + Gy(t). (38D) 
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APPENDIX E 


SAMPLE COMPUTER PROGRAM FOR OBSERVER/CONTROLLER DESIGN 


A computer program is given on the following pages which 
can be used to design an observer/controller in the manner 
prescribed by the algorithm of Section IID. It 1s assumed 


that the plant state and output equations are given by 


x(t) 


y Ce) 


Ax(t) + Bu(t) 


Cx(t) (1E) 
and that the observer state and estimate equations are given 
by 

z(t) = (A-GC)z(t) + Gy(t) + Bu(t). (2B) 


The control law is given by 
u(t) = Kz(t) (3E) 


and the performance measure to be minimized is given by 


r= f t¥pxTceyqx(t) + eT (t)Re(t) Jat. (4B) 
O 


The following symbols are explained in the order of their 
appearance in the input section of the program: 

Bed eee cmos PCO hen Ee tOm Ene Minimizing Search 
moter ne DD IRG I - 

EPS - The migtimum step length allowed. 

) = The minimum value that thé€tsum of the cigenvalues 


(trace) of the closed-loop system is allowed to take on. 


(AS 





a 





ALPH - The weighting multiplier which is used to assign 
a cost to any violation of the trace constraint. 

hie eine morderso themiderix Ain (LE). 

JY - The maximum number of searches the minimizing search 
routine is allowed to perform. 

NR - The time multiplier exponent r in the performance 
measure (4E). 

Man ic mEowedinecnsicn OF Coin (15)) For full-order 
observer/controller design M is set to zero. 

Aeoo ine Macri x A an) (VE); 

Q - The matrix Q in (4E). 

Bee hich matrix Bin (1E). 

GG - A vector whose elements are given by the initial 


estimate of G and K in the form 


So: The R Matrix in (4E). 
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